Effects of atomic diffraction on the Collective Atomic Recoil Laser 



OO 

on 



M. G. Moore and P. Meystre 
Optical Sciences Center and Department of Physics 
University of Arizona, Tucson, Arizona 85721 
(March 10, 1998) 

We formulate a wave atom optics theory of the Collective Atomic Recoil Laser, where the atomic 
center-of-mass motion is treated quantum mechanically. By comparing the predictions of this theory 
with those of the ray atom optics theory, which treats the center-of-mass atomic motion classically, 
we show that for the case of a far off-resonant pump laser the ray optics model fails to predict the 
linear response of the CARL when the temperature is of the order of the recoil temperature or less. 
This is due to the fact that in this temperature regime one can no longer ignore the effects of matter- 
wave diffraction on the atomic center-of-mass motion. 

PACS numbers: 42.55-f,42.50.Vk,03.75.-b 
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I. INTRODUCTION 

The Collective Atomic Recoil Laser, or CARL, is the 
atomic equivalent of the Free Electron Laser [Q. Devel- 
oped theoretically by Bonifacio et al [Q-|5| , the CARL de- 
vice has three main components: (1) the active medium, 
which consists of a gas of two-level atoms, (2) a strong 
pump laser which drives the two-level atomic transition, 
and (3) a ring cavity which supports an electromagnetic 
mode (the probe) counterpropagating with respect to the 
pump. Under suitable conditions, the operation of the 
CARL results in the generation of a coherent light field 
(the probe) due to the following mechanism. First, a 
weak probe field is initiated by noise, either optical in 
the form of spontaneously emitted light, or atomic in 
the form of density fluctuations in the atomic gas which 
backscatters the pump. Once initiated, the probe com- 
bines with the pump field to form a weak standing wave 
which acts as a periodic optical potential (light shift). 
The center-of-mass motion of the atoms on this poten- 
tial results in a bunching (modulation) of their density, 
very much like the combined effects of the wiggler and the 
light field leads to electron bunching in the free-electron 
laser. This bunching process is then seen by the pump 
laser as the appearance of a polarization grating in the ac- 
tive medium, which results in stimulated backscattering 
into the probe field. The resulting increase in the probe 
strength further increases the magnitude of the stand- 
ing wave field, resulting in more bunching followed by an 
increase in stimulated backscattering, etc. This positive 
feedback mechanism results in an exponential growth of 
both the probe intensity and the atomic bunching. This 
leads to the perhaps surprising result that the presence 
of the ring cavity turns the ordinarily stable system of an 
atomic gas driven by a strong pump laser into an unstable 
system. 

The operation of the CARL was verified experimen- 
tally by Bigelow et al ||], using a hot atomic cell. Re- 
lated experiments by Courtois et al using cold cesium 
atoms, and by Lippi et al Q using hot sodium atoms mea- 
sured the recoil induced small-signal probe gain, which 



was interpreted in terms of coherent scattering off an in- 
duced polarization grating. However, these experiments 
lacked a probe feedback mechanism, which is necessary 
to see the long time scale instability which characterizes 
the CARL. 

The CARL theory developed by Bonifacio et al consid- 
ers the atoms either as classical point particles moving in 
the optical potential generated by the light fields, or, in 
a "hybrid" version, as particles whose center-of-mass is 
labeled by their classical position, but with quantum fluc- 
tuations about that position included. From an atom op- 
tics point of view, such theories can be described as "ray 
atom optics" treatments of the atomic field, in analogy 
with the ordinary ray optics treatment of electromagnetic 
fields. 

Like ordinary ray optics, the ray atom optics descrip- 
tion of CARL is expected to be valid provided that the 
characteristic wavelength of the matter-wave field re- 
mains much smaller than the characteristic length scale 
of any atom-optical element in the system. The char- 
acteristic wavelength of the atomic field is its de Broglie 
wavelength, determined by the atomic mass and the tem- 
perature T of the atomic gas. The central atom-optical 
element of the CARL is the periodic optical potential, 
which acts as a diffraction grating for the atoms, and has 
the characteristic length scale of half the optical wave- 
length. Hence the classical "ray atom optics" descrip- 
tion is intuitively expected to be valid provided that the 
temperature is high enough that the thermal de Broglie 
wavelength is much smaller than the optical wavelength. 
This gives the condition T Tr, the recoil tempera- 
ture of the atoms, as the domain of ray atom optics. In 
particular, it is certainly expected to hold under the tem- 
perature conditions of the experiments performed so far. 

However, the spectacular recent progress witnessed by 
atomic cooling techniques makes it likely that CARL ex- 
periments using ultracold atomic samples can and will be 
performed in the future. In particular, subrecoil temper- 
atures can now be achieved almost routinely. The pur- 
pose of this paper is to extend the CARL theory to this 
"wave atom optics" regime In this regime matter- 
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wave diffraction is expected to play a dominant role in 
the CARL dynamics, and thus it becomes important to 
determine to what extent it counteracts the bunching 
process in the CARL. 

The wave optics theory of the CARL is similar to the 
analysis of atomic diffraction by standing waves |fC[ | , ex- 
cept that the electromagnetic field is now treated as a dy- 
namical variable. It is also similiar to the theory of recoil 
induced resonances [jllf] , which describes the stimulated 
scattering of light off a standing wave induced polariza- 
tion grating, but the absence of a feedback mechanism 
for the probe feedback in that case means that it lacks 
the instability necessary for lasing. 

In this paper we focus on the case of a far off-resonant 
pump laser, thus permitting us to neglect the excited 
state population and therefore to ignore the effects of 
spontaneous emission (except as a hypothetical source of 
noise for probe initialization). We further concentrate 
on the linear regime, where both the probe field and the 
atomic bunching are considered as infinitesimal quanti- 
ties, since it is this regime that determines whether or 
not the exponential instability occurs. Finally, we re- 
strict our analysis to atomic densities low enough that 
collisions between atoms may be ignored, and neglect 
the transverse motion of the atoms, which in the absence 
of collisions is decoupled from the longitudinal degree of 
freedom along which bunching occurs. 

We note at the outset that our theory is semiclassi- 
cal in that it treats the electromagnetic field classically. 
While this approximation can not fully describe the sta- 
tistical properties of the CARL output, it is sufficient 
to describe the small-signal gain of the system, provided 
that one makes the implicit assumption that small fluctu- 
ations will trigger it, an approach familiar from conven- 
tional laser theory and nonlinear optics. We also empha- 
size that it is not inconsistent to treat the matter waves 
quantum mechanically while treating the light classically, 
since the limits under which a quantum description is re- 
quired are independent. For light, this limit is usually 
associated with weak intensities, while for matter waves 
it is normally a low temperature limit. 

This rest of this paper is organized as follows: Section 
II briefly reviews the ray atom optics model of the CARL, 
establishing the notation and setting the stage for a com- 
parison of its predictions with those of the wave atom op- 
tics theory, which is introduced in section III. Section IV 
discusses the collective instability leading to CARL oper- 
ation, compares the ray atom optics and the wave atom 
optics predictions, and determines the domain of validity 
of the former theory. Finally section IV is a summary 
and outlook. 



al. [|]-f| It begins with the classical iV-particle Hamilto- 
nian 

JY 

H N = Y / H 1 (z j , Pj ), (1) 

where Zj and pj are the classical position and momentum 
of the jth atom, obeying the canonical equations of mo- 
tion d,Zj/dt — dH^/dpj and dpj/dt — —dH^jdzj. The 
single-particle Hamiltonian Hi is given explicitly by 

P TlLO 

Hi{zj,Pj) = ^ + -^a zj + ih [g x a\e- lk ^a^ 

+ g 2 a* 2 e- ik ^a^-c.c], (2) 

where m is the atomic mass, ojq is the natural frequency 
of the atomic transition being driven by the pump and 
probe lasers, and g\ is the atom-probe electric dipole cou- 
pling constant. It is given by gi = [i\\ck\j (2He V)] 1 ^ 2 , 
where p\ is the projection of the atomic dipole moment 
along the probe polarization, ki is the probe wavcnum- 
ber, and V is the quantization volume. The atom-pump 
coupling constant g 2 is defined analogously to gx, but 
depending on fi 2 , the projection of the atomic dipole 
moment along the pump polarization, and k 2 the pump 
wavenumber. The normal variables a\ and a 2 describe 
the probe and pump laser fields, respectively. They obey 
Maxwell's equation 

d N 

—on = -iujiai + gi 2J e~ lkiZi (T-j , (3) 
j-i 

where Ui is the natural frequency of the probe (i = 1) 
or the pump (i = 2) field. Note that these equations are 
also valid for quantized electromagnetic fields, provided 
that at are interpreted as annihilation operators, but we 
describe the light fields classically in this paper. 

The variables er_j and a Z j are the expectation val- 
ues of the quantum mechanical Pauli pseudo-spin oper- 
ators which describe the internal state of the jth atom. 
They obey the familiar optical Bloch equations, appro- 
priately modified to include the center-of-mass motion of 
the atoms and with spontaneous emission neglected [], 

4f-j = -iu <r-j + [g*a 1 e ikl ^ + g^a 2 e tk2Z '} a zj , (4) 
and 

^-a zj = -2 [ gi a* ie -^ + g 2 a* s e- lk ^] H + c.c. (5) 



II. RAY ATOM OPTICS MODEL 

The Ray Atom Optics (RAO) model of the CARL has 
been developed and extensively studied by Bonifacio et 



Spontaneous emission is neglected in anticipation of the 
future approximation that the pump lasers are far-off reso- 
nant, and therefore the excited state population may be safely 
neglected. 
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It is convenient to introduce slowly varying variables 
via the transformations ai — a' 1 e'^ wt , a 2 — a' 2 e~ lut , and 
(T-j = a'_.e~ l ^ t ~ k2Z i\ where uj is the pump frequency 
shifted by the frequency pulling contribution due to the 
induced atomic polarization. The exact value of uj will 
be derived in a self-consistent manner shortly in a way 
similar to the approach used in conventional laser theory. 
These new variables obey the equations of motion 
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dt 



rn 
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+ c.c. 
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-a; = i(w - wi)oi + 51 V e-^-^^V. 



dt 



3=1 



—a' 2 = i(u> - uj 2 )a' 2 



N 
3=1 



(6) 



(7) 



(8) 



(9) 



= -2 

+ c.c. 
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a /* e -»(fci-fc3)^ + 



/ * 
52 a 2 



(10) 



and 



d / k 2 
-a_ J = t { u - U0 --p J )a_ j 



9i 



/ e i(*i-fei)»j +5 * a ^ 



(11) 



In the case where the lasers are tuned far off reso- 
nance, and the atoms are initially in the ground state, 
the excited state population remains small and can be 
neglected. This is equivalent to describing the atoms as 
classical Lorentz atoms, and is accomplished by setting 
a z j = — 1 in Eq. (|ll|). Assuming further that the detun- 
ing uj — luq is much larger than any other frequency in 
Eq. (|i"l|), allows one to adiabatically eliminate cr'_j with 



{uj - lu ) 



9i 



la[e i{kl - k2)z > + 



92 a 2 



(12) 



where we have in addition neglected the Doppler shift 
k2Pj/m compared to uj — ujq. This leads to the reduced 
set of equations 



dt 



d 2hko 

—Pi = —i-, r 

at {uj — uj ) 



rn 

,* „ J * J „i2k Zi 



[g 1 g 2 a 2 a x e 3 - c.c. 



(13) 



(14) 



dt 



N\gi\ 2 
{uj - uj ) 



■ 9 2 9i 
(uj - UJ ) 



N 

a ' 2 e 

3=1 



- UJ\ 



-i 2 koZj 



(15) 



and 



dt 



N\g 2 \ 2 

(uj - UJq) 



UJ 2 



N 



(16) 



where we have introduced fco = {k\ — k 2 )/2. 

We now introduce the undepleted pump approxima- 
tion, valid in the linear regime where a[ remains small. 
This is achieved by dropping the term proportional to a[ 
in Eq. @. This yields 



—a 9 

dt 2 



N\g 2 \ 2 

(uj - UJq) 



Ul 2 



'2' 



(17) 



which has the steady state solution a' 2 (t) — a 2 (0) pro- 
vided that the frequency pulling condition 



N\g 2 \ 2 

{uj - UJq) 



-uj 2 = 0. 



(18) 



is satisfied. Note that this equation has two solutions, 
but we must choose the branch which gives the result 
uj = uj 2 when N = 0. This leads to the solution 



1 

U= 2 



oj + uj 2 ± \J(u) 2 -uj a ) 2 + AN\g 2 \< 



(19) 



where the plus sign must be taken for positive detunings 
{u>2 > ujq) and the minus sign for negative detunings 
{uj 2 < uJo)- Expanding this relation to lowest order in 
{uj 2 — ujq)^ 1 gives the expected result 



■ uj 2 + 



{uj 2 - lu ) ' 



(20) 



To proceed analytically past this point, it is convenient 
to introduce the dimensionless variables 9j = 2k^Zj 1 Pj — 
Pj/hko, A = gig2a 2 {0)ai/[uj r (uj — ujq)] and t = Auj r t, 
where the recoil frequency ui r is given by 

uj t = fikl/2m. (21) 

These variables obey the equations of motion 

d 



dr 



-6, = P. 



±P- 
dr J 



-iAe lHj + c.c, 



(22) 



(23) 



and 



3 



d 

rf7 



1 N 

A = iAA - ia— 

n 



-iB-i 



(24) 



3=1 



where we have introduced the dimensionless control pa- 
rameters 



and 



A = 5i/4w r , (25) 



- An<7i| 2 |<?2| 2 M0)|78^( w - Wo ) 2 , (26) 



where <5i = w — u>i — iV|<?i| 2 / (oj — u>o). We note that both 
A and a are real numbers, and furthermore that a > 0. 

We seek solutions of these equations which are per- 
turbations about the case A = 0. Thus we make the 
substitutions 



9 j = 6 j (0) + P j (0)T + 50 jl 



(27) 



and 



and 



dr 



-A 



AA-aY^B(k) 



(34) 



where N(k) is the number of atoms in the velocity group 
with momentum hkoP(k) and we have assumed that 



JV 

i=i 



-i20j(O) 



0, 



(35) 



an assumption that requires that N(k) ^> 1. Note that 
this formulation implies a discretization of the initial mo- 
mentum distribution, and furthermore assumes that the 
atomic positions in each velocity group are initially ran- 
domly distributed along the CARL cavity. Fluctuations 
in the initial distributions can of course readily be in- 
cluded into the initial conditions of the perturbation vari- 
ables. 



(28) 



where 9j(0) is randomly taken from a uniform distribu- 
tion, and Pj (0) is randomly taken from the initial mo- 
mentum distribution. The new variables 89 j and 8Pj 
give the perturbations on the atomic center-of-mass mo- 
tion due to a nonzero A(0). We introduce finally the lin- 
earized velocity group bunching parameter and its "con- 
jugate" momentum according to 



B(k) 



and 



N 



3=1 



N 



(l-^.) e -^i(0)+P 3 (0)r) ) (2Q) 



3=1 

+ P{k)B(k). 



-i(e j (0) + P J (0)r) 



We note that 



k 



B(k) 



-i2koz\ 



(30) 



(31) 



and the amplitude of (|31j) is a measure of the degree of 
bunching of the atomic gas. A magnitude of zero indi- 
cates no bunching, while a magnitude of one indicates 
maximum bunching. This leads to the equations 



— n{k) = i 



—B(k) = -in(k), 



P 2 (k)B(k) - 2P(k)n(k) - 



(32) 



N(k) 
N 



A 



(33) 



III. WAVE ATOM OPTICS MODEL 

In order to quantize the center-of-mass motion of a gas 
of Bosonic atoms, one may either utilize first quantiza- 
tion, and replace the variables Zj and pj in the TV-particle 
Hamiltonian ([!]) with operators satisfying the canonical 
commutation relations [zj,Pji] = ihSjj>, or equivalently 
we can second-quantize the single particle Hamiltonian 
(0), introducing creation and annihilation operators for 
excited and ground state atoms of a given center-of-mass 
momentum. It is this second method which we will adopt 
in deriving the Wave Atom Optics (WAO) model. In the 
absence of collisions, the second-quantized Hamiltonian 
is simply 



H = J2H(k), 



(36) 



where H(k) is given by 

h 2 k 2 



H(k) 



-&Uk)c g (k) 



h 2 k 2 



hu) I cl(k)c e (k) 



2m " gy "' HV " / ' V 2m 
+ ifi [gia\cl{k + k\)c e {k) + g2a* 2 c\(k + k 2 )c e (k) 
- H.c] , (37) 

where the field operator c g {k) annihilates a ground state 
atom of momentum hk, and c e (k) annihilates an excited 
atom of momentum hk. We assume that the atoms in 
the sample are bosonic, so that these operators obey the 
commutation relations 



[c g (k),cl(k')] = [c e (k),ct(k')]=S kkl 



(38) 



all other commutators being equal to zero. 

With the atomic polarization now expressed in terms of 
field operators, Maxwell's equations (||) for the classical 
laser fields become 



4 



dt ai 



-iwidi + gi)(cl(k + ki)c e (k)). 



(39) 



Hence, all that is required to determine the field evolu- 
tion are the expectation value of bilinear combinations 
of atomic creation and annihilation operators. The evo- 
lution of these expectation values is easily obtained by 
introducing the "single-particle" atomic density opera- 
tors 



p gg (k,k') = cl(k')c g (k), 



(40) 



— — p(k, k) — —i{n 2 — k' 2 )p(k, k!) 
dr 

+ ^A*[ P (k,k' + 1)-p(k-1,k')} 
-lA[p(K + l, K ')-p(K, K '-l)}, (45) 



and 



—A = iAA - ia^2p(n, k + 1), 



(46) 



p eg (k, k') = [p ge (k', fe)]t = cj (k')c e (k), (41) 

and 

p ee {k,k') = ct{k')c e {k). (42) 

Note that e.g. the expectation value of the diagonal op- 
erator (p g g(k 7 k)) gives the mean number of ground state 
atoms with momentum hk. The expectation values of 
these operators obey the equations of motion 



—Pjj,(k, k') 



~([H,p jf (k,k>)}) 



(43) 



where pjj>(k,k') = (pjj>(k,k')). The full form of these 
equations is given in the Appendix. The important point 
is that they are depend only on pjj>(k,k'), hence they 
form a closed set of equations which describe the response 
of the atomic field to the driving laser fields. We note 
that had we included collisions in our model, this would 
no longer be the case. 

Introducing in analogy to the ray optics description 



the rotating variables a\ 



, a 2 



a 9 e 



p eg (k,k') = p' eg (k — k 2 ,k')e~ lult , neglecting the excited 
state population, and solving adiabatically for p' eg (k,k') 
yields 



where the parameters A and a are given by Eqs. (|25| ) 
and (prjj), respectively. 

As in Sec. II, we seek a solution which is a perturbation 
about the case A — 0. From Eq. (pOf), the unperturbed 
solution is readily found to be 



p(k, k', t) = p(k, k', 0)e 



' nv~ i ( K2 ~ K ' 2 ) T 



(47) 



We consider specifically an atomic sample initially in 
thermal equilibrium, so that Eq. (071) becomes 



P(k,k ,t) = — — d KjK ,, 



(48) 



where N(k), the number of atoms with initial wavenum- 
ber 2kgK, is given by a thermal distribution function. We 
introduce the perturbation variables Sp(n, k') according 
to 



K ') = — + Sp(K, K'). 



(49) 



and anc ^ observe that Maxwell's equation (f46j), which becomes 



—A = iAA — ia Sp(n, k + 1), 



(50) 



Peg(k,k') 



(w - UJ ) 

gla 2 Pgg(k,k')\ 



[9i a 'iPgg(k + 2fc o, k') 



(44) 



Substituting Eq. (^4|) into Maxwell's equation (|39| ) for 
the pump and making once more the undepleted pump 
approximation leads to the solution a' 2 {t) = a 2 {0) pro- 
vided that lo is given by Eq. (|lS|). We then substitute 
Eq. ( |44| ) into the equation of motion for p gg {k,k'), and 
introduce the dimensionless wavenumber k — kj{k\ — k 2 ) 
and the mean density p(n, k') = p gg (k, k') /N , in addition 
to the dimensionless variables already defined in the ray 
atom optics model. We arrive at the wave optics equa- 
tions of motion 



2 These are single-particle operators in the sense of many- 
body theory, since they only involve the annihilation of an 
atom in a given state and its creation in some other state. 



together with the linearized equation 



—5p{n, k + 1) = i{2n + l)8p(n, n+l) 
dr 

[*(«+!) -*(«)] 

2N V ' 

form a closed set of equations which underlies the dy- 
namics of the CARL in the linear regime of wave atom 
optics. 



IV. COLLECTIVE INSTABILITY 

The most important feature of the CARL is the ap- 
pearance of a collective instability, which gives rise to 
exponential gain under appropriate parameter settings. 
This instability is characterized by an imaginary fre- 
quency component in the spectrum of the probe field 
A(t). As has been demonstrated in Ref. 0, one needs 
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not solve the complete set of equations derived in the pre- 
vious sections in order to determine the necessary con- 
ditions for the collective instability. Instead, by taking 
the Laplace transform of these equations one can derive a 
"characteristic equation" which allows one to determine 
whether exponential gain occurs, and if so what the ex- 
ponential growth rate is. 

For the Ray Atom Optics model, the Laplace transform 
of Eq. Q yields 



Ms) 
where R(s) is given by 



R( s y 



R(s) 



s — iA — ia 



f(k)dk 
{s + i2kf 



(52) 



(53) 



In obtaining this result we have taken the continuum 
limit and assumed that B(k) and TI(k) vanish at r = 0. 
Here f(k) is simply the normalized thermal distribution 
function for the dimensionless momentum P(k) = 2k. 
The roots of R(s) give the characteristic exponents of the 
CARL. Stability requires that all roots be purely imag- 
inary. When the collective instability occurs, however, 
there will be one root with a positive real part. This real 
part is the RAO exponential growth rate Tr. This result 
is identical to that first obtained by Bonifacio et al. 

The Wave Atom Optics model, which includes the ef- 
fects of atomic diffraction, yields the Laplace transform 



A w (s) 



MO) 
W(s)' 



(54) 



W(s) is given by 



W(s) 



iA — ia 



f(k)dk 



(s-i(2k-l))(s-i(2k+l)) 

(55) 



where we have again taken the continuum limit and as- 
sumed that 5p(k,k + 1) vanishes at r = 0. If a root of 
W(s) with a positive real part exists, that real part is the 
WAO exponential growth rate T^. 

We see by comparing Eqs. (52) and (54) that the effect 



of atomic diffraction is to lift the degeneracy of the sin- 
gularity under the integral. This expression also leads us 
immediately to the conclusion that if the width of the mo- 
mentum distribution f(k) is large compared to 2k, then 
the singularity will appear as essentially degenerate, and 
the effects of matter waves diffraction will be negligible, 
Thus the RAO and WAO models should agree for large 
enough temperatures. 



A. Finite temperatures 

In the absence of quantum degeneracies, the ther- 
mal momentum distribution is given by the Maxwell- 
Boltzmann distribution 



\/7T 



(56) 



where 1 — Tr/T and Tr = huj r /kB is the recoil tem- 
perature, Kb being the Boltzmann constant. By sub- 
stituting Eq. ( |56"| ) into Eq. ( |53| ) and using the Fourier 
convolution theorem we find that the RAO exponential 
growth rate is determined by the equation 



iA 



ia 



pe 



-p 2 /4(3 2 -ps 



dp = Q, 



(57) 



which can be integrated to give the transcendental equa- 
tion 

s-iA- i2a(3 2 + i2y/^a(3 3 eP 2s2 erk(f3s) = 0. (58) 

In contrast, substituting Eq. (|5^) into Eq. (]55|) and again 
using the convolution theorem we find that the WAO 
exponential growth rate Tw is determined by equation 



■ - iA - 



,-p 2 /40 2 -p S 



sin(p) = 0, 



(59) 



By examining Eq. (|59j) we see that in the case fi <C 1 
we are justified in expanding sin(p) to lowest order in p. 
This exactly reproduces Eq. (^), thus showing that the 
WAO and RAO descriptions make indistinguishable pre- 
dictions about the exponential growth rate in the limit 
T 3> Tr. However, for temperatures comparable to or 
less than the recoil temperature, we will see that the 
RAO theory fails to correctly predict the behavior of the 
CARL in the linear regime. Physically, this is due to 
the fact that it does not account for the effects of atomic 
diffraction, which tends to counteract the bunching pro- 
cess. Finally, we note that upon integration, Eq. (|5T 
becomes the transcendental equation 



3 l2/32s erfc[/3(s + %)] 



erfc[/?(s — i)\ 



0. 



(60) 



In the next subsection we will examine in more detail the 
precise manner in which diffraction interferes with the 
bunching process for the special case of a zero tempera- 
ture atomic gas. But before turning to this extreme situ- 
ation, we present numerical results comparing RAO and 
WAO models at non-zero temperature, as determined by 
solving Eqs. (||) and @. 

Figures l(b-d) compare Tr with Tw at a = 10 for the 
three different temperature regimes, T — Tr, T = IOTr, 
and T = IOOTr respectively. Figures 2(b-d) shows the 
same comparison for a = 10 -1 . While we see that the 
behavior of and Tw depends strongly on a (recall 
that a is proportional to both the pump intensity and 
the atomic density), the discrepancies between the two 
models as a function of temperature are very similar. 
At T = Tr there are significant differences between the 
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predictions of the RAO and WAO models, but these dif- 
ferences become minimal at T — IOTr, and insignificant 
at T = IOOTr. We also observe that the differences are 
more pronounced for lower values of a, meaning that 
at lower densities and/or pump intensities, the quantum 
mechanical behavior becomes more apparent. The reason 
for this is that at high intensities the bunching process, 
driven by the probe field, dominates, while at low inten- 
sities the anti-bunching effects of atomic diffraction play 
a larger role. 



B. The T — limit 

For a typical atom, the recoil temperature is of the or- 
der of microkelvins, e.g. for sodium we have Tr = 2ApK. 
However, recent advances in cooling techniques have 
led to measured temperatures as low as the picokelvin 
regime. At these extreme temperatures the condition 
T <C Tr is satisfied, i.e. we are effectively in the T — > 
limit . In this section we study the T = case in detail in 
order to gain further insight into the exact role of matter 
wave diffraction in the CARL system. 

For the RAO model, we have a single velocity group at 
k = 0. Thus by differentiating Eq. ( p2\ j with respect to r 
and using Eq. (B3J), we see that the bunching parameter 
B = -B(O) = (cxp(— i2k z)) obeys the equation of motion 

^B = -A, (61) 

where we have taken P(0) = and N(0) = N to indicate 
that all atoms are initially at rest. 

In the WAO description, setting N(k)/N = 5 K fl in Eq. 
(|l|) shows that two variables are coupled to the probe 
field, 8p(— 1,0), and 8p(0,l). They describe the recoil 
of atoms initially at rest as a result of their interaction 
with the light fields. We proceed then by introducing 
the new variable B = Sp(—1, 0) + Sp(0, 1), which has the 
same physical meaning as in the RAO model, namely 
B = (exp(— i2koz)). But in contrast to that case, the 
time evolution of B is now governed by the equation of 
motion 

^B = -B-A. (62) 

This result shows that in contrast to the predictions 
of classical mechanics, where the bunching parameter 
B has dynamics similar to a free particle driven by 
the probe field A, quantum mechanically B behaves 
as a simple harmonic oscillator of frequency 4w r (in 
original time units), and subject to that same driv- 
ing force. In the linear regime, B is assumed to be a 
small perturbation about its initial value of zero, and 
the forces resulting from a non-zero probe field A tend 
to cause B to increase. But this mechanism is op- 
posed by the "restoring force" due to atomic diffraction. 
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FIG. 1. Comparison of the exponential growth rate as a 
function of pump-probe detuning A between the RAO (solid 
line) and the WAO (dashed line) models, for the case a = 10. 
(a) shows the results for T — (see Sec. Ill.b), (b) shows the 
case T = Tr, (c) shows T = 10T«, and (d) shows T = 100T fl . 
We see that the ray atom optics model gives the correct result 
only in the limit T » Tr. 
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In addition to opposing any increase in the magnitude 
of B, the diffraction term also modifies its phase, which 
may upset any phase relation between A and B which 
might be required for the collective instability to occur. 

The RAO model only makes accurate predictions at 
T = in the limit lo t . — ► 0. Therefore, if we were to 
increase the mass of the atoms, thus decreasing w r , the 
behavior at T = would become more and more classi- 
cal. This is because heavier atoms suffer less diffraction 
than lighter atoms under the influence of the light fields. 
We also note that the correspondence principle states 
that quantum mechanics should agree with classical me- 
chanics in the limit Ti — > 0, which would also cause io r to 
tend to zero. These considerations can also be derived 
from the statement that the RAO model is valid when 
T 3> Tr, if we note that as u> r — ► the recoil temperature 
also goes to zero. 

In both the RAO and WAO models, the probe field A 
obeys the equation 



—A = i(AA-aB). 
uT 



(63) 



From this equation, together with Eq. ( pl\ ) we find that 
the solutions are exponentials with exponents given by 
the roots of the cubic equation 



iAs 2 



0. 



(64) 



This is exactly the "cold-beam" cubic equation of Bonifa- 
cio et al j|] . However with the inclusion of atomic diffrac- 
tion effects, we now see that the correct "cold-beam" cu- 
bic equation, derived from Eqs. ( |63| ) and (|6^), is 



s 3 - iAs 2 + s - i(a + A) = 0. 



(65) 



These equations can also be derived from the Laplace 
transform method of Sec. Ill, with the substitution 
f(k) = 5(k), indicating a zero temperature momentum 
distribution. 

From these cubic equations it is possible to determine 
the point of transition between the stable and the un- 
stable regimes of the CARL. For the RAO model the 
collective instability occurs provided that the threshold 
condition 



a > 



4A 2 



(66) 



is satisfied, and above threshold the exponential growth 
rate is given by 



FIG. 2. Figure 2 is identical to Figure 1, except we have 
now taken a = 10" 1 . Since a gives the stength of the bunch- 
ing process, when it is small the effects of atomic diffracton 
play a larger role, leading to stronger discrepancies between 
the predictions of wave atom optics and ray atom optics. 
However, we see that the RAO limit, given by T > Tr, is 
independent of a. 



r R = 



_ \/3 /a 



) 1/3 |(l + VC) 2 / 3 -(l-\/C) 2 / 3 |, (67) 



2 V4- 



where C = l-4A 3 /27a. For the WAO theory the thresh- 
old condition is 



a> 27 



(3 + A 2 ) 3 / 2 -9A + A 3 



(68) 
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and above threshold the exponential growth rate is given 
by 



\/3 /a\ 
~2~ V4 J 



1/3 



(1 + /D) 2 



2\2 



(l-^) 2 + ^(l-A 2 ) 2 



27a 2 

1/3 



(1-A 2 ) 



where 
D 



X _4A 

3a 



A 2 

1 
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2\ 2 



27a 2 



(l-A 2 ) 



1/3 



(69) 



(70) 



In Figure 3(a) we examine the CARL operating 
regime, defined as the region in parameter space where 
the exponential instability occurs, at T = as it 
would be if Ray Atom Optics were valid. We con- 
trast this with Figure 3(b) which shows the actual 
CARL operating regime at T = 0, as calculated us- 
ing Wave Atom Optics. From this figure we see that 
the operating regime of the CARL is drastically re- 
duced at low pump intensities and/or atomic densi- 
ties when the effects of atomic diffraction are included. 

20 




FIG. 3. The CARL operating regime (shaded region) as 
predicted by the RAO model (a), and the actual operating 
regime (b) , as given by the WAO model. 



see that this second threshold may even be above A = 
for low intensities and/or densities. In fact, the threshold 
crosses A = at precisely a = 2/3V3- 

Figure 2(a) shows that in the limit of weak pump inten- 
sities and/or atomic densities the peak gain for the WAO 
model tends to A = 1, while that of the RAO model is 
at A = 0, This result can actually be understood quite 
simply: The atomic center-of-mass dispersion curve tells 
us that the absorption of a pump photon and the emis- 
sion of a probe photon by an atom initially at rest cre- 
ates an energy defect of 4w r due to atomic recoil. This 
defect can be compensated by a detuning between the 
pump and probe, which in dimcnsionless units occurs at 
A = 1. Therefore, the fact that T\y is a sharply peaked 
function around A = 1 is simply an expression of energy- 
momentum conservation. If we are to take the Ray Atom 
Optics model seriously at T = 0, then we must concede 
that we are in the limit where u r — > 0, therefore, energy- 
momentum conservation would predict the maximum of 
Tr to occur at A = 0. In other words, in that limit the 
center-of-mass dispersion curve is flat over the range of 
a few photon momenta. 



V. DISCUSSION AND OUTLOOK 

The main result of this paper is that at low tempera- 
tures the behavior of the CARL is strongly influenced by 
matter-wave diffraction, which tends to counteract the 
atomic bunching and reduces the instability range of the 
system. In this temperature range, the CARL presents 
an experimentally realizable example of dynamically cou- 
pled Schrodinger and Maxwell fields. The present theory 
quantizes the matter wave, but not the electromagnetic 
field. It will be of considerable interest to extend it to 
regimes where both fields need to be quantized. An anal- 
ysis of the density regime where quantum degeneracy be- 
comes important will also be a fascinating extension, in 
particular when two-body collisions are included. This 
study will allow one to investigate to which extent a Bose- 
Einstein condensate can be manipulated and modified in 
a far off-resonant CARL configuration. An intriguing 
possibility would be to generate in this fashion a cou- 
pled laser- "atom laser" system. The study of the coher- 
ence properties of this system will be the object of future 
investigations. Finally a comparison between bosonic 
and fermionic CARL systems in the quantum degenerate 
regime should also be considered. 



ACKNOWLEDGMENTS 



Figure 1(a) compares Yr with for the case a = 10 
at T = 0, and Fig. 2(a) shows the same comparison 
for a = 10 _1 . We see that atomic diffraction leads to 
the appearance of a second threshold below which the 
collective instability does not occur. From Fig. 2(a) we 



We have benefited from discussions with R. Bonifacio 
and L. De Salvo, who brought the CARL system to our 
attention. This work is supported in part by the U.S. 
Office of Naval Research Contract No. 14-91-J1205, by 
the National Science Foundation Grant PHY95-07639, by 



9 



the U.S. Army Research Office and by the Joint Services 
Optics Program. 



APPENDIX: HEISENBERG EQUATIONS OF 
MOTION FOR DENSITY OPERATORS 



The full equations of motion for the expectation values 
of the density operators are 



V)p 99 (k,k>) 

+gialp eg (k - ki,k') + g 2 a* 2 p eg (k - k 2 , k') 
+9iaip ge (k, k' - fci) + g 2 a 2 p eg (k, k' -k 2 ), 



(Al) 



d_ 

dt 



Peg (k, k ) 



±(k?-k> 2 )+u; 
Ira 



Peg (k, k ) 



+gla\[p ee {k,k' - h) - p gg (k + k\,k')] 
+g 2 a 2 [p ee (k, k' - k 2 ) - p gg {k + k 2 ,k% (A2) 



and 



4>ee(M') - -^(k 2 - k' 2 )p ee (k,k') 

at 2m 

-g\a\p eg {k,k' + fci) - g 2 a 2 p eg (k,k' + k 2 ) 
-glaip ge {k + ki,k') - g 2 a 2 p ge (k + k 2 ,k'). (A3) 
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